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Abstract 

In this paper we introduce the well-balanced Levy driven Ornstein-Uhlenbeck pro- 
cess as a moving average process of the form Xt = J exp(— A|t — u\)dLu- In contrast 
to Levy driven Ornstein-Uhlenbeck processes the well-balanced form possesses con- 
tinuous sample paths and an autocorrelation function which is decreasing not purely 
exponential but of the order A|m| exp(— A|u|). Furthermore, depending on the size of A 
it allows both for positive and negative correlation of increments. We indicate how the 
well-balanced Ornstein-Uhlenbeck process might be used as mean or volatility process 
in stochastic volatility models. 
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1 Introduction 



Recently moving average processes have attained much attention, both from the theoretical 
and application side, since they provide a large class of processes, only partly belonging to the 
class of semimartingales and allowing to model correlation structures including long-range 
dependence. The theoretical foundations of treating moving average processes with driving 
Levy processes have been provided in Rajput and Rosinski (1989) and recently the question 
under which conditions these type of processes are semimartingales has been considered 
in Basse and Pedersen (2009). A special case of Levy driven moving average processes are 
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fractional Levy motions (cf. Benassi et.al (2004) and Marquardt (2006)), where the kernel 
function of the fractional Brownian motion is taken, leading to the same correlation structure 
as fractional Brownian motion. Bender et.al (2010) derived conditions on the driving Levy 
process and the exponent of the kernel function under which the fractional Levy motion is a 
semimartingale. It turns out that this can only be the case in the long memory setting and 
then the process is of finite variation. Barndorff-Nielsen and Schmiegel (2009) developed the 
idea of moving average processes further by introducing a stochastic volatility component 
leading to Brownian semi-stationary processes, which are a very promising class of processes 
for modelling turbulence. Furthermore, these processes have also been applied to electricity 
modelling (cf. Barndorff-Nielsen et.al (2010)). However, we can also view the well-known 
Ornstein-Uhlenbeck process as moving average process, which due to its simple structure is 
very popular for modelling mean reverting data (e.g. Barndorff-Nielsen and Shephard (2001), 
Kliippelberg et.al (2009)). 

Motivated by this we introduce an exponential kernel exp(— A|t — -j), A > on the whole 
real line leading to the well-balanced Ornstein-Uhlenbeck process. We show that this process 
is well defined with only assuming a logarithmic moment on the driving Levy process. The 
process possesses infinitely divisible marginal distributions and is stationary. In contrast 
to Levy driven Ornstein-Uhlenbeck processes it possesses continuous sample paths of finite 
variation and therefore it is a semimartingale with respect to any filtration it is adapted 
to. Furthermore, the autocorrelation function is decreasing more slowly than the one of 
the Ornstein-Uhlenbeck process with same A, namely it is of the order X\h\ exp(— A|/i|). In 
addition the range of the first-order autocorrelation of the increments is (—0.5, 1) in contrast 
to (—0.5, 0) for the Ornstein-Uhlenbeck process. Positive values are often associated to long 
range dependence, but with the well-balanced Ornstein-Uhlenbeck process we see that this 
is not true. 

Hence the well-balanced Ornstein-Uhlenbeck process might serve as a promising mean pro- 
cess in financial models, e.g. as additive component in stochastic volatility models, since it 
possesses the following desirable properties: 

• the decay of the autocorrelation function is of the order X\h\ exp{—X\h\), 

• the autocorrelation between increments can be positive and negative, depending on A, 

• it is a semimartingale, 

• it has an infinitely divisible distribution. 

Furthermore, as in Barndorff-Nielsen and Shephard (2001) the well-balanced Ornstein- 
Uhlenbeck process might be used as volatility process in stochastic volatility models. Since 
the relationship of cumulant transforms between price and volatility process only differs by 
a constant from the Ornstein-Uhlenbeck model, modelling of the marginals stay the same 
as in Barndorff-Nielsen and Shephard (2001). However, the different correlation structure of 
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the well-balanced Ornstein-Uhlenbeck process is inherited by the integrated volatility and 
the squared price increments. 

In addition to the well-balanced Ornstein-Uhlenbeck process with the kernel given above we 
also introduce the process with the corresponding difference kernel exp(— A|t— -1)— exp(— A|- 1), 
motivated by the form of the kernel of fractional Brownian motion. This process, in contrast 
to the previous one, is not stationary, but it possesses stationary increments and starts in 
zero. Furthermore, the distribution of the squared increments of both processes are obviously 
equal and the autocorrelation function has the same decay. 

Let us give a brief outline on how the paper is organized: in Section 2 we introduce the nota- 
tion and define the processes, in Section 3 we show that both processes are semimartingales 
and derive the structure of their characteristics. In Section 4 we provide the moments and 
correlation structure of the processes. In Section 5 we give a brief empirical example to high 
frequency data. In Section 6 we indicate how the well-balanced Ornstein-Uhlenbeck process 
might be used as volatility process. 



2 Definition of the well-balanced Ornstein-Uhlenbeck 
process 

As driving process we consider a Levy process L given by the characteristic function 
E{exp{iuLt)) = exp(tip{u)) with 

iIj{u) = iuj — h / ( exp{iux) — 1 — iuxl\x\<i ) i'{dx), 

2 J -oo ^ ^ 

where the Levy measure z/ satisfies the integrability condition 1 A vi^dx) < oo. 

In the following we give conditions on a kernel function /(■,■) : x M — )■ such that 
processes of the form 

/oo 
f{t,s)dLs, t>0 
-oo 

exist. Here L denotes the two-sided version of the Levy process which is defined in the 
straight forward way by taking two independent copies L^^^ and L^'^^ and defining 

lJ^^ if t > 
-I'll if t < 0. 

Here and in the following we deal with stochastic integrals on the real line as well as on 
the positive half line. Integrals on M are meant in the sense of Rajput and Rosinski (1989), 
i.e. we associate an independently scattered random measure A with the two-sided Levy 
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process L. For details we refer the reader to Sato (2004) (in particular Theorem 3.2) who 
even treats the more general case of additive processes in law on [0, oo). The extension to 
M is straightforward. A is defined on the 5-i'm.g of bounded Borel measurable sets in M and 
the integral j^g{s) dAg is introduced in a canonical way for deterministic step functions g. A 
function / is then called integrable if there exists a sequence {gn)nm of step functions such 
that 



dn f a-s with respect to the Lebesgue measure 
lim„_^oo fAdni^) ^-^s exists for every A G 



If a function / is integrable, we write J^fdLs = limn^oo J^gn{s) dAg. From time to time 
we will switch between this integral and the classical Ito integral, namely in the case 
Jjg l[o,t]/(s) dLg = /g /(s) dLg for / G Ci,. Both integrals coincide for predictable integrands 



of the type /(s) 



-[0,s) 



. The general case follows by a standard argument using dominated 



convergence. Before we specify the function /(t, s) let us first briefly look at the setting of 
a general kernel. Rewriting the criteria of Rajput and Rosinski (1989) [Theorem 2.7] for the 
existence of the integral we obtain: the stochastic integral f dLg is well defined if for t > 



oo POO 



oo J — oo 



x/(t, s)|^ A 1 v(dx) ds < OO 

/oo 
a'^f{t,sfds < oo 

/OO /'OO 
/(t, 5)^7+ / x(l|^/(t,^)|<i - l|^.|<i) z/((ix) 
oo J —oo 



ds < oo 



(cf. in this context Basse and Pedersen (2009) equations (2.1)-(2.3)). Then the characteristic 
function is given by 



E{exp{iuZt)) = exp ip{uf{t, s)) ds 
and Zt is infinitely divisible with characteristic triplet (7^, (t|, z/j) 

fit, s)(^ + j x{l\^f{^t,s)\<i - l|x|<i) v{dx)j ds 



If 



-00 
00 



aV(t, sfds 



Vf[A) = (z/x A)|(a;,s) x/(t,s) G A\{0}}, A e B 

where A denotes Lebesgue measure. Furthermore for ui,U2, - ■ ■ ,Um G M and — 00 < ti < 
t2 < ■ ■ ■ < tm < 00 we obtain 

m „ m 

E(^expC^iUjZt^)^ =exp(^ tp{^Ujf{tj,s))ds 
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If we now consider kernels of the form f[t — s) the resulting process Z is stationary since for 
every h > 

m „ m 

E(^expC^iujZt^+h)^ = exp M i'{^Ujf{tj + h - s)) ds^ 

/m 
i=i 

m 

i=i 

In particular it possesses stationary increments. If we consider kernels of the form f{t — s) — 
f{s) we have Zq = a.s. and stationary increments where the increments have the same 
distribution as the increments of the process generated by the kernel f{t — s). 

If f{t, .) G L^(M) and the second moment of L exists and the first one vanishes, we denote 
E{L\) = V, then Zt also exists in the L^-sense with isometry 

EZ^ = \\f{t,.)\\l,V, 

as shown in Marquardt (2006) [Prop. 2.1]. Now we can come back to our special cases and 
assume A > 0. For the well-balanced Ornstein-Uhlenbeck process the kernel is 

exp(— A|t — s\) = exp ^ — A( max(t — s, 0) + max(— (t — s), 0)) j . 

From this reformulation we can see why we call the process well-balanced Ornstein-Uhlenbeck 
process, namely 

/oo /»^ /"OO 

exp{—X\t — s\)dLs= / exp{—X{t — s))dLs+ / exp{—X{s — t)) dLs 
-oo J ~oo Jt 

which is analogous to the well-balanced fractional Levy motion (cf. Marquardt (2006) [Defi- 
nition 3.1]). Following the terminology of Samorodnitsky and Taqqu (1994) [Example 3.6.4] 
in the stable case it is the sum of an Ornstein-Uhlenbeck process and a reverse (or fully 
anticipating) Ornstein-Uhlenbeck process. The initial distribution of X is given by 

/•oo 1*0 i*oo 

Xo= e~^l"lrfL, = / e^'dL,+ / e'^' dL^. 



oo 



As for the fractional kernel we can construct processes with stationary increments starting 
from zero, which for the well-balanced Ornstein-Uhlenbeck process leads to 

/oo 
exp(-A|t-s|) -exp(-A|s|)rfL,. 
-oo 

Now we can provide the characteristic triplet of the process X. 
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Lemma 2.1. The well-balanced Ornstein-Uhlenbeck process 

/CO 
exp(— A|t — s\) dLg 
-oo 

is well-defined and infinitely divisible with characteristic triplet {'^x,<^X','^x] 
Ix = t7 + T / '^{dx) - / u{dx] 



1 
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<rx = J- 

Px{A) = (z/ X A)|(x,s) xexp(-A|t- s|) G A\ {0}|, A e 13 
if and only if X > and J A log |x|z/((ix) < oo. 

Proof. The result follows by straight forward calculations from the general formulae. □ 

Here we see that in contrast to fractional Levy motions the well-balanced Ornstein-Uhlenbeck 
process is well-defined with only imposing the condition of a logarithmic moment of the 
driving Levy process. 

3 Semimartingale Property and Characteristics 

Since the processes X and Y differ only by a random variable which does not depend on 
t > 0, in the following we only treat X. However, the results remain valid for Y. 

We will show that {Xt)t>o is a process of finite variation on compacts and hence a semi- 
martingale with respect to any filtration it is adapted to. In order to do this we introduce 
the following decomposition 

/oo fO ft POO 

e-^l*-^l dLs = e"^* / e^^ dL^ + e"^* / e^^ dLs + e^* / e"^^ dL^ 
-oo J —oo Jo Jt 

and write the last term as 

l*00 POO pt 

e^* / e"^" dLs = e^* / e"^" dL^ - e^* / e"^" dL^. 
Jt Jo Jo 

For reference purposes we write the above representation of {Xt)t>o in a short form 

Xt = e-^'G + e^'H + e~^'lt-e^'Jt, (1) 
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using the following notation: 



t nt 

As Jr „„J T . / „— As 



I^ ■= I e^' dLs and Jt := / e~^' dL^ 
Jo Jo 

and 

(•0 POO 

G := I e^' dL, and H := / e'^' dL,. 





From ([T]) above and 1.4.36 in Jacod and Shiryaev (2003) we obtain: 

AXt = A(e-^*/t - e^'Jt) = e-^\e^'ALt) - e^\e-^'ALt) = 

for every t > 0. This means the process is continuous, though the driving term exhibits 
jumps. The following theorem refines this result. 

Theorem 3.1. The process {Xt)t>o is of finite variation on compacts and hence it is a 
semimartingale with respect to any filtration it is adapted to. Furthermore the process is 
locally Lipschitz continuous. 

Proof. First we use the integration- by-parts formula on If and —Jf and obtain 

t j-t 

e^' dLs = - LsXe^' ds + L^e^* 
Jo 

respective 

ft 

e-^' dLs = - Ls\e~^' ds - Lte~^K 
Jo 

Putting these together we have 

/O roc rt nt 

e^' dLs + (e^* - 1) / e"^* dL^ - e"^* / L,Ae^' ds - e^* / L.Ae"^" ds 
oo Jo Jo Jo 

PS PS 

+ Xe-^' / LrXe^'' dr - e-^'L^Xe^' - Xe^' / L.Ae"^" dr - e^U^Ae"^" J> ds 
Jo Jo 



- GXe-^' + HXe^' - 2L,A + Ae"^" / LrXe^' dr - Xe^' / L.Ae"^" dr } ds 
1 Jo Jq 

X{ - Ge-^' + He^' - I^e'^' - J^e^' )> ds 



A 



-Ar 



dLr — e 



-As 



dLr ) ds 



This means we end up with a pathwise integration of functions which are a.s. locally bounded 
since they are cadlag. In particular the process is locally Lipschitz and of finite variation on 
compacts. □ 

Note that this is different to the classical Ornstein-Uhlenbeck process which inherits the 
jump property from the driving process. 

In Basse and Pedersen (2009) and Bender et.al (2010) the authors treat the case of other 
kernel functions. However, their conditions on the Levy process are more restrictive. Note 
that by Jacod and Shiryaev (2003) Proposition 1.4.24 the process (Xt)f>o, since it does not 
have jumps, is even a special semimartingale with respect to every filtration it is adapted 
to. For the remainder of the paper we fix the filtration F which is obtained by defining first 
= {J^})t>Q via J^} := a{J-'^,G,H). Which is completed and made right continuous in the 
usual way to obtain F. 

Since the process is continuous and of finite variation on compacts we obtain the following 
corollary by the general representation of semimartingales (cf. Jacod and Shiryaev (2003) 
Theorem II. 2. 34). 

Corollary 3.2. The semimartingale characteristics of the process X are {Xt — Xo,0,0). 

A different approach leads to another perspective on the well-balanced Ornstein-Uhlenbeck 
process. The process X can be represented as an ARMA process in continuous time, i.e. 
a CARMA process (cf. Brockwell and Lindner (2009)). In general such a process, possibly 
complex valued, is given by 

Xt = h'Rt 

where b E and the C^-valued stochastic process R is given as the solution of the SDE 



dR, 



1 





1 





■ 


\ 














1 
























Rtdt + 















1 













— ap„i 






) 







dLt 



where ai, Cp are complex- valued coefficients and L is a driving Levy process. The p x p- 
matrix is usually denoted by A. In our case we have 



p 



2, A 



1 
A2 



-2A 




The SDE for the well-balanced Levy driven Ornstein-Uhlenbeck process is 



dR 



(1) 



dt + 
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dRr = X'^Rt dt + dLt 

with the exphcit solution 

Rt = e^'Ro + dK. 

Here e"^* is not meant componentwise but in the sense of bounded operators, i.e. as a con- 
vergent series of matrices. By choosing the initial conditions appropriately, i.e. 

we get our process X as a stationary CARMA(2,0) process. For criteria for CARMA processes 
to be stationary compare Brockwell and Lindner (2009) [Theorem 3.3]. 

Summarizing we can see that though integrating with respect to a quite general Levy process 
the special very regular form of the kernel leads to a semimartingale of bounded variation. 
Hence the well-balanced Ornstein-Uhlenbeck process might serve as mean process in the 
framework of semimartingale models, e.g. stochastic volatility models in finance. 



4 Moments and Correlation Structure 

In this section we will analyze the correlation structure of the well-balanced Ornstein- 
Uhlenbeck process. We will see that though the process is closely related to the stationary 
version of an Ornstein-Uhlenbeck process the two-sided kernel leads to a different behaviour 
in the autocorrelation function, namely to a decay of the order A|/i| exp(— A|/i|), and to a 
bigger range of possible values as in classical Ornstein-Uhlenbeck process, including positive 
ones, in the first order autocorrelation of increments. 

Proposition 4.1. Let Xf = J exp(— A|t — u\)dLu and assume that the driving Levy process 
possesses a finite second moment. We denote the variance of L{1) by V and the first moment 
by fi, then we obtain the following characteristic quantities for X and h > 



EXt 
var{Xt) 



2/i 
T 



Cov{Xt+h,X,) = Vhe-^'^ + je-^'' 
Corr{Xt+h,Xt) = A/ze"^' + e-^^ 

Proof. From the general form of the characteristic function, we can calculate the second 
moment of = f f(t,s)dLs, provided that L possesses a second moment and both f(t, .) 
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and f{t, are integrable. We obtain 



EZt = j f{t,s)ds(^ + j xl\j:\>iv{dx)^ 

EZl = j f{t,sfds (^a'^ + j x^u{dx)^ + (^j f{t,s)ds^ ^ j ^^l^l>i^(^^)) • 

In the following we denote o"^ + j x^u^dx) = V and 7+ J xl|a;|>ii^((ix) = fi. Using this together 
with the independent increment property of L, we obtain for Xf = J exp(— A|t — u\)dLu and 
s < t 



Hence 



EX, 

var{Xt) 



2/i 
T 

A A2 
V 



A 

EiXt-X^y = ^(l-e-^(*-^)-A(t-s)e-"(*-^)). 
A 



Cov{Xt,Xs) = ^{EX^ - E{Xt - Xsf + EX^) - EXtEXs 

= V{t - s)e-^(*-^) + ^e-"(*-^) 

A 

Corr{Xt,Xs) = Xit - s)e~''^'~'^ + e~^^'~'\ 

□ 



Comparing this to the well known quantities of a stationary Ornstein-Uhlenbeck process U 
we see, while the mean and the variance only differ by a multiple of two, the auto-covariance 
and autocorrelation function have an extra term leading to a slower decay, when the same 
A is considered. This might be an interesting feature for modelling data, especially coming 
from finance, where a pure exponential decay often seem too fast to match the empirical 
autocorrelation properly. 

If we define the well-balanced Ornstein-Uhlenbeck process with time scaled by A, i.e. X, = 
J exp(— A|t — u\)dLxu this has the same effect as in Barndorff-Nielsen and Shephard (2001), 
that the marginal distribution and hence the moments are independent of A. 

Also the correlation between increments might be of interest for modelling purposes and 
follows by direct calculations from the proposition above. 
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Corollary 4.2. Assume the same conditions on L as in the previous proposition, then we 
obtain 

1 1 1 - exp(A) + Aexp(A) 



Corr{Xk+i - Xk, Xi - Xq) = exp{-Xk) 

+Xk exp{—\k) 



2 2 1 -exp(-A) - Aexp(-A) 

1 1 1 - exp(A) + Aexp(-A) 

2 ^ 21-exp(-A)-Aexp(-A) 



and as a special case the first- order autocorrelation 

'1 + A 1 1 + A-exp(A) + A2exp(-A) 



Corr{X2 - Xi, Xi - Xq) = exp(-A) 



2 1 — exp(— A) — A exp(— A) 



Note that in contrast to the classical Ornstein-Uhlenbeck process whose autocorrelation 
function of increments Corr{Uk+i — Uk,Ui — Uq) = exp{—Xk){^ + ^jz^^~X)) always 
negative in the range between -0.5 and 0, we can have positive and negative values, in 
the range from -0.5 to 1 depending on A for the well-balanced Ornstein-Uhlenbeck process. 
Looking for example at the first-order autocorrelation Corr{X2 — Xi, Xi — Xq) it is positive 
for A < 1.25643 and negative for bigger values of A. This provides much more flexibility for 
modelling, e.g. we can obtain values for the first-order autocorrelation which is often linked 
to long-range dependence. Assuming that Bj^ denotes a fractional Brownian motion with 
Hurst parameter He (0, 1), then by Kettani and Gubner (2006) 



^Ch = 22^-1 - 1 



where X/^ = Bf^ — B^ ^ and X^ = - . X/^. Hence we can see that as the first-order 
autocorrelation of the well-balanced Ornstein-Uhlenbeck process Ch G (—0.5, 1) and Ch > 
for H > 0.5. 

For some applications it might of course be more realistic not to have a stationary process, 
but a process with stationary increments like Levy processes. In the context of well-balanced 
Ornstein-Uhlenbeck processes we can construct processes with the same correlation structure 
of increments and the same paths regularity by considering the associated difference kernel. 

Proposition 4.3. LetYt = exp(— A|t — s|) — exp(— A|s|)(iLs and assume that the driving 
Levy process possesses a finite second moment. We denote the variance by V and the first 
moment by fi, then we obtain the following characteristic quantities for Y 

EYt = 

var{Yt) = Vte^^' + ^e-^' 

A 

^ f ^r\^^ 1 1 - exp(A) + Aexp(A) A 

Corr{Yk+i - n, Y, - Fq = exp{-Xk) - + f\{ , ^ \. > 

\2 2 1 — exp(— A) — Aexp(— A)^ 



\z, z i — expi^— a; — A expi^— /\j y 

, /I 1 1 -exp(A) + Aexp(-A) ^ 

+Xkexp(-Xk) - H ^t^-^t , \ ,\ 

^ V2 21 -exp(-A) - Aexp -A 
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Proof. The proof follows immediately by noting that Yt = Xt — Xq. 



□ 



Note that we can easily also construct a process which only possess this correlation structure 
for a specific lag and is zero for larger lags. For a kernel on a compact interval [0, a] we 
obtain the process Xt = J^*_^exp(— A(t — s))dLs which possesses the second moment EX^ = 
(1 — exp(— 2Aa))/(2A), assuming for simplicity that the first moment of -^^(1) is zero and the 
second 1. Furthermore for increments Xt — Xs we obtain E{Xt — XsY = (1 — exp(— 2Aa) — 
exp(-A(t-s))+exp(-A(2a+s-t)))/A iit-s < a and iit-s > a: E{Xt-X,f = EX^+EX^. 
This leads to Cov{Xt, Xg) = (exp(— A(t — s)) — exp(— A(2a + s — t)))/(2A) for t — s < a and 
otherwise. 



5 Application to high frequency data 

We apply the well-balanced Ornstein-Uhlenbeck process to an example of real data and 
show that the autocorrelation models the empirical autocorrelation quite well. Hence this 
indeed offers the possibility of adding the well-balanced Ornstein-Uhlenbeck process X as 
an empirically convincing mean process to a classical stochastic volatility model, i.e. for the 
log-price process Yt = Xt + /J asdWs + Jt, where a denotes a volatility process, W a Brownian 
motion and J a jump Levy process. 

First we consider one trading day of the SAP share, namely of 1st February 2006 9:00 am 
to 5:30 pm consisting of 5441 trades. 



autocorrelation function SAP SAP 




200 400 600 800 1000 200 400 600 800 1000 

lag frequency 
5441 trades 



The picture shows the empirical autocorrelation function as solid line, the dashed line is the 
fit with a classical Ornstein-Uhlenbeck process and the dotted line with the well-balanced 
Ornstein-Uhlenbeck process. We can see that the autocorrelation function both visually and 
by taking the residual sum of squares fits the data much better than the Ornstein-Uhlenbeck 
process, except for small lags. This might be interpreted as the effects of market microstruc- 
ture. Namely the two kinks in the empirical curve are at a lag of 75 and 150 respectively. 
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In this setting this correspond to a samphng frequency of 7 minutes and 14 minutes. Val- 
ues in this range are in the econometrics literature often seen as samphng frequencies from 
which market microstructure effects start to be negligible. We have also included the volatil- 
ity signature plot for the one trading day of 1st February 2006, i.e. the realized volatility 
Sr=i(^i/" ~ -^(j-i)/n)^ plotted with decreasing sampling frequency. It shows the presence 
of market microstructure when the realized volatility increases with increasing sampling fre- 
quency. For the SAP data we can see indeed this effect, which vanishes above roughly 150, 
but possesses an additional peak at about 80. If we start fitting the empirical data only for 
larger lags than 150, the values of A and the RSS for the Ornstein-Uhlenbeck process stay 
the same, whereas the RSS of the well-balanced Ornstein-Uhlenbeck processes decreases to 
4.39310-1 

Furthermore, we also examined different data sets, namely tick-by-tick data of the Daimler 
Chrysler share and the Siemens share of January 2005. Fitting the empirical autocorrelation 
with the autocorrelation function of the well-balanced Ornstein-Uhlenbeck process and the 
Ornstein-Uhlenbeck process up to a lag of 1000, the well-balanced Ornstein-Uhlenbeck did 
better on 17 of the 21 trading days for Daimler Chrysler and on 13 for Siemens. Hence for 
single days the well-balanced Ornstein-Uhlenbeck process provides a good alternative. For 
averaged data the situation is a bit different. When averaging of the trading days the well- 
balanced Ornstein-Uhlenbeck process performed better than the Ornstein-Uhlenbeck process 
for Daimler Chrysler and worse for Siemens. Especially for small lags the Ornstein-Uhlenbeck 
process seems to fit very well as averaging smoothes away the kinks in the empirical curves 
of the single days. 

6 Well-balanced Ornstein-Uhlenbeck Process as 
Volatility Process 

In this section we examine how the well-balanced Ornstein-Uhlenbeck process might serve 
as volatility process in a stochastic volatility model along the lines of the Non-Gaussian 
Ornstein-Uhlenbeck models of Barndorff-Nielsen and Shephard (2001). The well-balanced 
Ornstein-Uhlenbeck process leads to a model which is as tractable as the Ornstein-Uhlenbeck 
model but possesses some new features, namely that the different decay of the autocorrelation 
function is inherited by the decay of the autocorrelation of increments of integrated volatility 
and the autocorrelation of squared log-returns. 

As in Barndorff-Nielsen and Shephard (2001) we consider the stochastic volatility model of 
the type 

dYt = {a + l3Xt)dt + xl''^dWt. 
We assume a, /3 G M, ly denotes a Brownian motion and X the spot volatility process given 
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by 

/oo 
-oo 

with A > and L a two-sided Levy process. We assume that the law of L{1) is self- 
decomposable and restricted to the positive half line and the cumulant generating function 
is given by 



k{e) = \og{E{exp{~eL{l)))) = - 1 - exp{-ex)u{dx), 

Jo+ 

where z/ denotes the Levy measure and the first two moments of L(l) exist. We denote the 
first moment by fi and the variance by V. The time scale At in the definition of X ensures 
that the marginal law is independent of A and the support of p that X is non-negative as 
required for a volatility process. 

Similarly as for the characteristic function in Section 2 also the cumulant transform k of X 
can be expressed in terms of the kernel and the cumulant transform of the underlying Levy 
process, namely 

/■oo roo 

k{9) = / k{9e-\^'-'\)ds = 2 / k{9e-'')du. 



From this we obtain that the n-th cumulant m„ of the underlying Levy process is related 
to the n-th cumulant m„ of X by m„, = |m„. Furthermore, if u possesses a Lebesgue 
density g, also i^g-|.| possesses a Lebesgue density g{y) = 2 g{xy)dx or respectively g{y) = 
^{~9{y) ~y9'{y))- This implies that the results differ only by 2 and 0.5 from the formulae of 
the Ornstein-Uhlenbeck processes respectively, hence modelling of the marginals is the same. 

Next by direct calculation we can provide an explicit formula for the integrated volatility 

Xudu = j{2L,t + Xo" - X+ - (Xr - 

where Xt = e~^^^~'^^dLxu + G~^^^~^^dL\u = X^ + X^. Hence we can see that we have 
the additional terms ^(-^^a* — X^ +X^) compared to the formula of the Ornstein-Uhlenbeck 
process which is provided in Barndorff-Nielsen and Shephard (2001). 

Finally we would like to note that the correlation structure of X is inherited by the integrated 
volatility and the squared returns. By the general formulae provided in Barndorff-Nielsen 
and Shephard (2001) we obtain for n,s > 1 setting a = /3 = 

Cov I / Xudu, / Xudu I = VR{As) 

[n~l)A J{n+s-l)A 



Corr{(YnA - 'i^(n-l)A)^ {y{n+s)A - y(n+s-l)A)^) 

where for t > 



R{As) 



6f(A) + 2A^f 



r{t) = Xte-^' + e-^' 
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f{t) = [ [ r(x)rfxdn = ^(Ate"^* + 2At + 3e^^*-3) 
Jo Jo ^ 
R{As) = f(A(s + l))-2f(As)+f(A(s-l)) 

= ^e-^'^'liXAs + 3)(e"^^ + e^^ - 2) + AA(e-^^ - e^^)]. 
A 
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